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We consider the causal structure of generalized uncharged McVittie spacetimes with increasing 
central mass m{t) and positive Hubble factor H{t). Under physically reasonable conditions, namely, 
a big bang singularity in the past, a positive cosmological constant and an upper limit to the central 
mass, we prove that the patch of the spacetime described by the cosmological time and areal radius 
coordinates is always geodesically incomplete, which implies the presence of event horizons in the 
spacetime. We also show that, depending on the asymptotic behavior of the m and p[ functions, 
the generalized McVittie spacetime can have a single black hole, a black-hole/white-hole pair or, 
differently from classic fixed-mass McVittie, a single white hole. A simple criterion is given to 


distinguish the different causal structures. 
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I. INTRODUCTION 

Fully time-dependent solutions of general relativity 
may behave radically differently from their stationary 
counterparts. Familiar properties assumed for most so¬ 
lutions, many of which require exact time independence 
and asymptotic flatness at heart, do not hold on the more 
general dynamical scenarios [1-3]. This breakdown of 
stationary properties on non-stationary spacetimes jeop¬ 
ardizes the interpretation and our very understanding of 
generic features of familiar spacetimes that we usually 
take for granted. 

Moreover, the search for quantum gravity, as well as 
the resolution of the most pressing theoretical difficul¬ 
ties faced by cosmology, such as the cosmological con¬ 
stant problem and its many related issues, has spawned 
a rapidly expanding collection of modifications of gen¬ 
eral relativity, many of which can be interpreted in the 
Einstein frame as additional fields on top of the classical 
Einstein-Hilbert action. These theories have been found 
to yield several physically interesting non-vacuum solu¬ 
tions of the Einstein equations, many of which are by 
construction non-static. Even stationary solutions have 
been raising debates on such established results as the no¬ 
hair theorem [4], prompting for a revisit of such results 
and a renewed interest in checking its ranges of validity, 
as well as the precision of its statements. 

In the study of the causal structure of non-vacuum 
solutions of general relativity, many concessions need to 
be made when stationarity is fully abandoned. Global 
properties have to be forfeit in favor of local ones, and 
inferences on the asymptotic structure need to take the 
bulk behavior into account. 
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One such example of a class of exact solutions to a 
modified gravity problem whose properties can be ac¬ 
cessed via analytical means is the McVittie metric [5], 
which describes a black hole in an asymptotically FLRW 
universe. It was recently found [6] that it is a solution to 
the incompressible limit of fc-essence minimally coupled 
to general relativity, also known as the cuscuton field [7]. 
Some generalizations of the McVittie metric have been 
proposed in the literature, such as a charged central ob¬ 
ject [8] and a time-dependent central mass [2]. The lat¬ 
ter has become known as the generalized McVittie metric 
(hereby referred to as gMcVittie for brevity) and has also 
been shown to be an exact solution of Einstein’s equa¬ 
tions for a scalar field coupled to gravity [9], specifically 
a particular case of the Horndeski theory [10] which can 
also be interpreted as an imperfect fluid in which a radial 
heat flow accounts for the increase in the central object’s 
mass [11[. It is related to the original fixed-mass McVit¬ 
tie metric via a disformal transformation, in which case 
the transformed action belongs to the more general 
class of scalar fields [12] . 

The richness of the original McVittie solution’s causal 
structure has been explored in previous works. It was 
shown that this metric actually describes a black hole in 
the appropriate limits [13, 14], and also that the history 
of the free functions present in the metric is responsible 
for the asymptotic behavior [15] . In the generalized case 
the situation becomes much more complex due to the 
addition of one more almost arbitrary time dependency 
coming from the mass function, and the analysis needs to 
be done without the simplifying properties of a constant 
central mass. 

In this work we determine whether the physically in¬ 
teresting limits of the gMcVittie metric are geodesically 
incomplete and examine the number and nature of the 
apparent and event horizons of this spacetime. We find 
that when the expansion tends asymptotically to de Sit¬ 
ter as time progresses, and when the central mass tends 
to a finite upper value, the gMcVittie spacetime is always 
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geodesically incomplete. Moreover, whereas the original 
McVittie metric had two distinct asymptotic structures 
depending on the expansion history, namely a black hole 
and a black-hole/white-hole pair, we find that the gen¬ 
eralized case presents a third possible outcome, in which 
only the white-hole horizon is present. 

The paper is organized as follows: in Sec. II we review 
the apparent horizon structure of the gMcVittie metric 
and show that under physically meaningful assumptions 
these spacetimes are geodesically incomplete with respect 
to ingoing null radial trajectories, which implies that the 
geometries have event horizons that can be associated 
with either a black hole or a white hole; in Sec. Ill we 
determine on which of these categories the event horizon 
falls, according to its dependence with the bulk history of 
the mass and expansion functions; we present our conclu¬ 
sions in Sec. IV. Throughout the paper, derivatives with 
respect to the t coordinate are denoted with an overhead 
dot. We use signature (—,-I-,-I-,-I-) and natural units 
with G = c = 1. 


II. STRUCTURAL ANALYSIS OF GMCVITTIE 
SPACETIMES 

A. General structure 


erties^. This corresponds to surfaces where 

0in0out = 0 , (5) 

where 0in/out the expansions of spherical congruences 
along radial null directions (the subscript “in” refers to 
the ingoing direction while “out” refers to the outgoing 
direction). Regions where 0in0out < 0 are called regular 
regions and regions where 0in0out > 0 are trapped re¬ 
gions if 0in/out < 0 or anti-trapped regions if 0in/out > 
0.2 

Setting ds2 = 0 for null geodesics in Eq. (1), it imme¬ 
diately follows that the relevant equations are 

( ^ ) = R{rH ±R)-^rMil-R) =0, (6) 

V dt / 

where, due to our coordinate choice, 0 < i? < 1. Since 
the turning points of 0in/out coincide with the roots of 
(®)in/out accreting 

black hole {M > 0) in an expanding universe {H > 0), 
only the ingoing null geodesics, which correspond to the 
minus sign, may have real roots in a finite (t,r) patch^. 
Therefore in the analysis of the apparent horizon struc¬ 
ture we only consider the minus sign. 

If we use the definition of R from Eq. (2) to substitute 
r as the main variable for which to solve Eq. (6), we can 
cast it as 


The uncharged gMcVittie metric, characterized by a 
mass function m{t) and a scale factor o(t), is given in 
areal radius and cosmological time coordinates by 


ds" = -R^t'^+i ^ - 


dr 


R 


H+M [--1 


rdf )■ -l-r^dfl^ ^ 


with 



(1) 

III 

1 

to 

(2) 


(3) 

m 

(4) 


We impose the condition that M < H at all times so 
that the causal structure of the singular surface r = 2m, 
namely the fact that it is spacelike everywhere and to the 
past of all causal curves, is retained from the standard 
McVittie metric [11]. The motivation for this somewhat 
restrictive condition is to keep the “McVittie big bang” 
a regular Cauchy surface from which to formulate a well 
behaved initial-value problem that allows us to use null 
geodesics to build the causal structure of the spacetime 
[14, 16[. 

Apparent horizons are defined as the surfaces in which 
congruences of null geodesics change their focusing prop- 


R^ -R? = -2m [{H - M)R-^M] . (7) 

This is a fourth-order algebraic equation, whose real so¬ 
lutions in the interval 0 < i? < 1 (r > 2m) correspond to 
the loci in which radial null geodesics have constant ra¬ 
dius in the ingoing direction. For continuous M and i7, 
a spacelike hypersurface containing an odd-multiplicity 
solution corresponds to the existence of a boundary be¬ 
tween an anti-trapped and a regular region, while even- 
multiplicity solutions separate regions which are either 
both anti-trapped or both regular. 

Since there is no cubic term in Eq. (7), all solutions 
must add up to zero, and there can be at most three pos¬ 
itive real solutions [22]. Moreover, Eq. (7) is helpful for 
the determination of the number of real roots admissible 


^ Here we are using the expression “apparent horizon” as a syn¬ 
onym of “trapping horizon” as coined in Ref. [17]. This coincides 
with the use of the term in Ref. [18], for example, but it is not 
equivalent to the definition given in Ref. [19], since the latter is 
not a quasi-local definition. 

^ This naming convention is found, for example, in Ref. [20]. Some 
authors prefer to use “normal” or “untrapped” instead of “regular” 
and “future trapped” (“past trapped”) instead of “trapped” (“anti¬ 
trapped”). 

^ This property is also observed in the Schwarzschild-de Sitter met¬ 
ric if one switches from isotropic coordinates to areal radius coor¬ 
dinates. The outgoing horizons are brought to finite coordinates 
via a subsequent transformation in the t coordinate, which the 
McVittie metric does not allow. 
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in the range 0 < i? < 1, which is the physically interest¬ 
ing one. If we think of the horizons as the intersection 
between the curve represented in the left-hand side with 
the straight line in the right-hand side, it becomes man¬ 
ifest that there can be at most two positive roots in the 
allowed range for R, as in the McVittie metric. For that 
reason, we adopt the same naming convention and de¬ 
note the inner apparent horizon as r_(t) and the outer 
horizon as r+(t). It can also be seen that, for certain 
histories of M and H, pairs of horizons may appear and 
disappear as the spacetime evolves, even if the restriction 
M < if is always respected [23]. Figure 1 illustrates this 
structure for one particular example of M and ii at a 
fixed time slice. 



Figure 1: Visualization of the horizons in the physically 
relevant patch of generalized McVittie as intersections 
between the curves associated with the left- and 
right-hand sides of Eq. (7). 


This horizon structure implies that, if we restrict 
our analysis to gMcVittie spacetimes that tend to non¬ 
extremal Schwarzschild-de Sitter as t —>■ oo, we are guar¬ 
anteed to have a regular region bounded by two ap¬ 
parent horizons, r_ and r_|_, for sufficiently large times. 
At large r, the geometry behaves asymptotically as 
a Friedmann-Lemaitre-Robertson-Walker (FLRW) met¬ 
ric, which means that all future-directed outgoing null 
geodesics extent to null infinity, since they always reach 
distances sufficiently large to be well approximated by 
outgoing geodesics in FLRW space. On the other hand, 
future-directed ingoing null geodesics that start in the 
normal region tend to fall towards r_, accumulating in 
its neighborhood as t —>■ oo at radius 

Too = lim r_{t) . (8) 

t—>-oo 

It is this limit surface {t —>■ oo,roo), mapped by ingoing 
null geodesics, that may constitute a black-hole horizon. 
If ingoing null geodesics departing from some initial event 
reach t —>■ oo in a finite affine parameter, it means that 
the patch of spacetime described by Eq. (1) is geodesi- 
cally incomplete. This issue has been thoroughly dis¬ 


cussed in the literature regarding the standard McVittie 
spacetime [13, 14[ and special cases of gMcVittie in which 
the time derivative of the mass function has support in a 
finite interval [11[, and in these cases it was shown that 
under certain conditions the surface (t — >■ oo,roo) consti¬ 
tutes a traversable limit surface. 


B. Geodesic incompleteness 


To write the geodesic equations in more physically 
meaningful terms, we use as a starting point the comov¬ 
ing flow in isotropic coordinates = —Rdt, which cor¬ 
responds to the Hubble flow in the asymptotically FLRW 
region. The expansion scalar 0(u) of tho timelike flow 
is given by 


©O/'i = 3 


H + M[ - 



= 3n{r,t), 


(9) 


so that H is defined as a generalization of the Hubble 
factor in the asymptotic analysis. We then rewrite the 
radial null trajectory in Eq. (6) as 


^^=-R{R-nr). (10) 

Using this result, the radial null geodesic equations may 
be cast as 


d^r 

dA2 

dH 

dW 


dt [R{R — Rr)] / dr \ 

R^{R-nr)^ ’ 


M + 


H - M 
2 



( 11 ) 


2m 

J.2 



2 


( 12 ) 


Since the surface (t —>■ oo, Too) is at a coordinate infin¬ 
ity, to assess the geodesic completeness of gMcVittie we 
need to check whether radial null geodesics starting from 
initial conditions at some event in the spacetime reach 
the surface (t —> oo,roo) after a finite affine parameter. 
If this is the case, this surface is a causal horizon and the 
spacetime can be analytically extended beyond it. Oth¬ 
erwise, the surface is at an infinite affine distance from 
any point in the spacetime, so it is not an event horizon, 
but another patch of null infinity (along with the patch 
defined by outgoing null geodesics). The late-time be¬ 
havior of ingoing null geodesics for large r allows us to 
establish the following theorem: 

Theorem II.1. The patch of gMcVittie solutions de¬ 
scribed by metric (1) in the {t,r) coordinates, with smooth 
m(t) and H(t) for all t > 0 and under the following hy¬ 
potheses: 


m{t) > 0, Vt > 0, 

(13a) 

Too = lim m{t) > 0, 

t—^OO 

(13b) 

Hq = lim H{t) > 0, 

t—yoo 

(13c) 
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3v^ ^ ^ ^ ’ 

(13d) 

M(t) > 0, Vt > 0, 

(13e) 

H{t)-M(t) > 0, VO 0, 

(13f) 


is null-geodesically incomplete. 

Each of the hypotheses grouped in Eqs. (13) of the 
theorem has a physical motivation. Eq. (13c) means 
that there is an expanding cosmological background; 
Eq. (13b) means that the mass is bounded; Eq. (13d) 
is interpreted as the spacetime having a non-extremal 
Schwarzschild-de Sitter limit as t —>■ oo; Eq. (13e) cor¬ 
responds to an accreting central object; and Eq. (13f) 
means that the spacetime presents a spacelike singu¬ 
larity in the past. Under these assumptions, Theorem 
II. 1 guarantees that all gMcVittie spacetimes are null 
geodesically incomplete and have a traversable surface at 
{t —>■ oo,roo). Note that Eq (13f) excludes the Sultana- 
Dyer spacetime [24], which corresponds to a gMcVittie 
space with H = M. As a particular case of this result, 
we include all standard McVittie spacetimes satisfying 
H{t) > 0, Hq > 0 and > mHo > 0, which general¬ 
izes"^ the results of [13]. 

Since ingoing null geodesics are uniquely determined 
by any two of Eqs. (10), (11) or (12), we choose to work 
with Eqs. (10) and (11). From Eq. (10), we note that in¬ 
going null geodesics r{t) always approach the inner hori¬ 
zon r-{t), since R — Hr > 0 above r_ (unless we are 
also above r+) and R — Hr < 0 below r_. This means 
that all ingoing radial null geodesics that contain events 
in the regular region at sufficiently large times tend to 
Too as t ^ oo, as well as those that are below r_ at suf¬ 
ficiently large times. Hypotheses (13b), (13c), (13d) and 
(13e) guarantee that there exists T such that we have 
the two apparent horizons r_ and r+ for all t > T. With 
those observations in mind we are ready to present the 
proof of Theorem II. 1, by computing the leading term of 
the expression for AA for an ingoing null geodesic taken 
near r^o in the limit as t —>■ oo. 


Proof. In order to study the asymptotic behavior, we can 
take t large enough so that we can linearize our expression 
near the t ^ oo limit value of the functions. For economy 
of notation, we denote as S the order of magnitude of our 
linear approximations: 


S 


max 



AHrao,Mroo, 



(14) 


We also define 7?oo = ^(wo,roo) and AH{t) = H{t) — 
Hq. Using this notation we can approximate R[m(t),r] 


at large times, up to o(6) terms, by 


-^oo H dynRooij^ 

d,.R^(^r Too) “f ^(^) 
TO — Too Too 


= Rco- 


oo-^n^ + ;:^(r-r^) + a(5), 


(15) 


and, since limt_,,oo M{t) = 0 because of Eq. (13b), we can 
also write 


rH{t,r) =rooHQ + r, 

+ (?' - roc)HQ + o{6). 


AH{t) + ^ (1 - Roo) 


(16) 


Recalling that R^o — r^Hg = 0 and ^ ^ ^°° , we 

obtain at first order in 6 


R-Hr = 


l-3Rg 

2i?oo 


r -roo \ TO - Too 
r oo Roo 


M 


AH+ —{I-Roo) 

-^oo 


(17) 


o(J). 


Near roo, Eq. (10) becomes at order 6 


dr 

dt 


= -Ro 




+ o{6) , 


(18) 


with 


i:U\ _ w - Too 

• oo-**^oo 


M 

AH P—{I-Roo) 


a = 




(19) 

( 20 ) 


In order to simplify our expressions, we also define z{t) = 
r{t) — roo analogously to Ref. [15] . In terms of z, Eq. (18) 
reads 


dz 

dt 


-uHqZ + Roo^{t) -f o{S), 


( 21 ) 


whose solution, with initial condition z(to) = Zq, is 


z{t) = i?ooe"“"^“‘ f + ^06”“"^“* + o(<5), 

Jto 

( 22 ) 

where Zq = zoe°‘^°*°. Now, working on Eq. (11) near 
r = Too or small z, we obtain at order 6 


dt [R{R - Hr)] 
R^{R-Hr)^ 



Ro 


ijt) 


+ o{6) , 


In [13] the additional assumption that H(t) < 0 was made. 


(23) 
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where 


m = 


H{t) + 

-^OO 


• (24) 


We obtain for z{X) 


d^z 

dW 


C(t) 


Bo 




dz 

dA 


+ o(S). (25) 


Integrating the leading term of Eq. (25) once, we find 


and changing the integration variable to t once 
more, we obtain at leading order 


KAA = Alloa 



Boo 

cxHqA 


i{t')dt' - aHot 


(32) 

Therefore, geodesic incompleteness is equivalent to 
the convergence of the expression 


dt. 



^ OO 

aA 




aHot 


dt , 


(33) 


dz 

-=ifexp 




[- 



(26) 


where iil is a constant and t should be understood as a 
function of z. 

Here we need to analyze the leading term of z(t), which 
depends on the behavior of ^(t). Three distinct regimes 
are possible: 


1. If ^{t) = o(e we can choose t large enough 

such that the solution reduces to 


z(t) = + o(e-“-f^«‘). 


(27) 


where the constant A codifies the dependence of 
z{t) on the initial conditions. 

2. If ^{t) > for large t, we have® 

1 


z(t) = — 
a 


m - 

Hoa 


o[e(t)]. (28) 


3. If C(i) = Cooe-°‘^°^ + o(e-“'f^«*), where C^o is a 
constant, we can solve the integral in Eq. (22) to 
find 


z{t) = i?ooCoote-“^“‘ + 0(e-“^“‘). (29) 


where we have used the fact that Boo = ^oo^o to 
further simplify the expression. Recalling that, in 
this case, ^{t) = then necessarily ^{t) = 

o(e““^“‘). This implies that J* ^{f)df < 


j* ^aHot ^-aHot go the leading term in 

the argument of the exponential in Eq. (33) is the 
linear one. Thus, Eq. (33) reads 


1= e-“^«*'df'+ o(e-“'f^“*), (34) 

Jto 

which is convergent. Therefore, the ingoing 
geodesics reach the limit surface in a finite param¬ 
eter time. 


Case 2: Here, ^(t) dominates the late-time behavior of 
the flow. The deduction works in the same way as 
for case 1 up to Eq. (26), but when we express the 
integrand as a function of t the expression behaves 
differently. Taking z = ^C(i) - + o(C) and 

dz = ^^{t)dt + o(^), we find for the leading term 


dz 

dA 


— K exp 


Bo 


a‘ 


'■m 




(35) 


We now work out each of these cases individually. 

Case 1: We may use Eq. (27) to change the inte¬ 
gration variable to t, which implies that dz = 
— -|- o(e““^°‘)] dt. Eq. (26) then 

reads at leading order 


- = /f exp 


Bo 




aHoA , 
Solving for A, we find 

Boo i 


i{t)dt 


KAX = - 


exp 


aHoA , 


^aHgtA 


^{t)dt 


dz. 


(30) 


(31) 


^ More detail about this approximation can be found in ap¬ 
pendix A. 


We then integrate again to obtain AA: 


pO 

K'A\= / 

J Zo 


^ CXD 

a 


^ — CxHnt 


i{t)dt, 


'to 


(36) 


so in this case geodesic incompleteness is equivalent 
to the convergence of the expression 


r= 

Jto 


dt)dt. 


(37) 


However, the integral (37) always converges when 
^ 0 as t —> OO, which we are assuming. 


Case 3: With ^ = CooC -I- o(e “^o*), we have z(t) 
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given by Eq. (29), so that 


dz = - [i?ooC'ooai?oie-“^“* + 0(6-“'^^“*)] dt. 

Thus, Eq. (26) takes the form 

dz ,, [ i?oo /■ m , ] 

[--e(^r 1 (38) 

= K', 

K'AX = -zo- (39) 

This finishes our proof that in all cases satisfying hy¬ 
potheses (13), the surface {t —>• oo,roo) is always reach¬ 
able by ingoing null radial geodesics in a finite affine pa¬ 
rameter. □ 


III. CAUSAL STRUCTURE 

The results of Sec. II state that accreting gMcVittie 
spacetimes which satisfy the hypotheses (13) can be ex¬ 
tended beyond the (t —> oo, Too) limit surface. The natu¬ 
ral question which arises in this situation is what kind of 
surface it is and what kind of region is hidden beyond it. 
In the standard McVittie scenario there has been a se¬ 
ries of articles discussing this issue [13-15], showing that 
in McVittie the {t oo,roo) limit surface can either be 
a black-hole event horizon or be composed of a black- 
hole horizon and a white-hole horizon. The distinction 
between these cases, as was shown in Ref. [15], is how 
fast the ingoing null geodesics approach Too in relation 
to the apparent horizon r_. If all geodesics approach 
the Too limit from above r_, then the limit surface is a 
black-hole event horizon. If, on the other hand, there 
are ingoing null geodesics that approach the limit sur¬ 
face from below r_(t) at the same time as ingoing null 
geodesics that approach the limit surface from above r_, 
then the limit surface contains a black-hole event horizon 
portion “above” r_ as well as a white-hole event horizon 
portion, “below” r_. These particular cases are depicted 
in Figs. 2, 3a and 3b, respectively®. 

This distinction comes from the fact that the Too sur¬ 
face in the Schwarzschild-de Sitter spacetime, which is 
the limit of gMcVittie as t —)■ oo given (13), is a zero 
of the expansion of outgoing null geodesics, ©out- Since 
©out > 0 in the patch described by the (f,r) coordi¬ 
nates, this means that behind the limit surface we have 
©out < 0- Then, the nature of the region behind the 
surface will depend on the sign of ©in. Since between r_ 


The model used to generate Figs. 2, 3a, 3b and 4 is H{t) = 
i/otanh“^ with Hg = 10“^, and m{t) = mg + 

Ata,nh{BHgt — 2.5). Fig. 2 uses mg = 5, A = 0.5 and B = 5; 
Fig. 3a uses mg = 8, T = 2, B = 5; Fig. 3b uses mg = 8, 
A = 0.025 and B = 1.5; and Fig. 4 uses mo = 10, A = 1 and 
B = 0.5. 


and we have ©;„ < 0 (it is a regular region), then if it 
is possible for timelike or null observers to travel to the 
limit surface while being always above r_ and reach it, 
then by continuity they will fall in a trapped region, as 
©out©in > 0 there. This is equivalent to saying that they 
have entered a black hole whose event horizon is the limit 
surface (or part of it). Conversely, below r_ we have an 
anti-trapped region with ©out > 0 and ©in > 0. If we 
can cross the limit surface while always traveling below 
r_, then by continuity we reach a regular region where 
©out < 0 and ©in > 0. This is equivalent to saying that 
we came out of a white-hole horizon. For more details in 
this argument, we refer the reader to Ref. [14] . 



Figure 2: gMcVittie spacetime where all ingoing null 
geodesics (in blue) reach the limit surface from above 
the horizon. 

In the following development, we will establish suffi¬ 
cient conditions over the functions defining each gMcVit¬ 
tie model, m{t) and H{t) that determine to which causal 
structures it corresponds. In this wider class of space- 
times, there is a new case, depicted in Fig. 4, in which all 
ingoing null curves approach the limit surface from below 
the apparent horizon r_. In this case, the limit surface 
is entirely a white-hole horizon. We can see in Figs. 2, 
3a, 3b and 4 that each case corresponds to a different 
end point of the r_ apparent horizon, when it reaches 
the Too limiting surface in the conformal diagrams. The 
causal structure is determined by the asymptotic behav¬ 
ior of the function ^{t), defined by Eq. (19) according to 
Theorem III.l. 

Theorem III.l. Let be a gMcVittie spacetime described 
by the metric (1) under the hypotheses (13) and, in ad¬ 
dition, assume that the function f,(t) defined in Eq. (19) 
tends to zero from positive (negative) values. Then, if 
there exists a > 0 such that 

FZ{to,t)= f (40) 

Jto 

converges as t ^ oo, then the gMcVittie spacetime 
presents a black hole and a white hole horizon. Anal¬ 
ogously, if there exists a > 0 such that 


F+i*o,t) 




■t 


(41) 
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Figure 3: gMcVittie spacetime where some ingoing null 
geodesics reach the Too limit below r_ (blue). Since 
later ingoing null curves emerge from the r = 2m 
singularity at the left, we see that there is a limit time 
that separates curves that reach Too from above from 
curves that reach there from below. 



Figure 4: gMcVittie spacetime where all ingoing null 
geodesics reach the limit surface from below the r_ 
apparent horizon (blue). 


diverges as t ^ oo, then the gMcVittie spacetime in ques¬ 
tion presents only a white hole (black hole) horizon. 


The results can be more easily visualized in Table L 
From Theorem IILl, one can see that the cases where 
gMcVittie contains only a white hole are restricted to 
cases where ^ 0+, that is, when it goes to zero from 

positive values. Symmetrically, the only black hole case 
correspond to ^ 0“. This is due to the fact that the 


Table I: Possible asymptotic structures of gMcVittie 
spacetimes: summary of the results of Sec. III. 



^{t) —)• 0 

m ^ 0+ 

g(affo ‘^)“^(u)du < OO 

black hole 
and 

white hole 

black hole 
and 

white hole 

JO® e(“^“+‘^^“.J(u)du —>■ 00 

black hole 
only 

white hole 
only 


sign of is the sign of the slope of r_ (t) for large t, 


r-{t) 


dr{R-Hr) 


r=r- (t) 


(42) 


since the denominator is positive at the r_ horizon. We 
are assuming that the denominator is not degenerate, 
that is, r_|_ and r_ do not coincide. We can understand 
intuitively that when its slope is negative, it is easier 
for null rays to reach the horizon from above the ap¬ 
parent horizon, which correspond to the regular region, 
which leads to a black hole. In the same manner, when 
the slope of r_ is positive, it is easier for null rays to 
traverse from below the apparent horizon, laying in the 
anti-trapped region, characterizing the limit of a white 
hole region. In both cases, if the absolute value of func¬ 
tion ({t) decreases fast enough, we have the case in which 
the limit surface correspond to a pair of white-hole/black- 
hole horizons, separated by a bifurcating two-sphere. The 
cases in which the limit surface has only one character 
are those in which the ^ function does not decrease faster 
than the exponential that modulates it in Eq. (41). 

With these remarks in mind, we can build models cor¬ 
responding to the four cases shown in Table I. First, we 
need to chose if the slope of r_ will be positive or neg¬ 
ative for large times, which by Eq. (42) means choosing 
the sign of ^ for large times. Inspecting Eq. (24), we can 
see that the term proportional to M is positive, by our 
initial assumptions. The sign of H can be either plus or 
minus, but physically realistic models usually correspond 
to H ^ 0. T^his means that v^e can tune the functions 
m(t) and H(t) in such a way that the leading term for 
large t is positive or negative. 

The second step is to choose how the leading term ap¬ 
proaches zero. If we choose m —toq or AH (depending on 
what was chosen to be leading before) that approaches 
zero sufficiently faster than exjp{—aHQt), then there will 
surely exist a such that FZ converges and we are in the 
case where a pair black hole/white hole is present. If 
they decrease sufficiently slower than exp(—aFfgt), then 
we fall in the case were diverges and there is only a 




















black hole or white hole, depending on the former choice 
of sign for ^ being negative or positive, respectively. Re¬ 
call that a and JIq depend only on the limit values of the 
gMcVittie functions, but the causal structure depends on 
how those limiting values are reached. Using the reason¬ 
ing above we have built the models depicted in Figs. 2, 
3a, 3b and 4. 

The proof of Theorem III.l is analogous to the proof 
of the main result in Ref. [15] and the case-specific de¬ 
tails can be found in Appendix C. Roughly speaking, the 
proofs for both cases are based in the fact that we can 
substitute r_(t) by r = r^a, which is a constant curve. 
Thus, we can evaluate the flow of ingoing null geodesics 
near r = r^o and determine if ingoing null curves cross 
or fail to cross it for sufflciently large t. This is done by 
defining two families of approximative curves, one that 
approximates the ingoing geodesics from above and an¬ 
other that does so from below. The existence or non¬ 
existence of approximating curves that cross the r = Too 
line for arbitrarily large times leads to the criteria shown 
in Table I. 


IV. CONCLUSIONS 

In this work we have shown that all generalized McVit- 
tie spacetimes with sufficiently slow accretion, an ex¬ 
panding cosmological background, a big bang singular¬ 
ity in the past and a non-extremal Schwarzschild-de 
Sitter spacetime as its limit as t —?> oo are geodesi- 
cally incomplete, having a traversable surface located at 
(t -)■ cxo,roo)- 

This result was proved by studying the leading term 
of the variation of the affine parameter of ingoing null 
geodesics as they approach the surface which accumu¬ 
lates near the inner apparent horizon r_(t). This was 
done by approximating the flow of the non-linear differ¬ 
ential equation that governs them for large times, taking 
care to keep the right leading terms at every step. We 
have shown that, under the aforementioned assumptions, 
the ingoing geodesics always reach the surface in a finite 
affine parameter, implying that the patch of generalized 
McVittie spacetime, described by the cosmological time 
t and the areal radius r in Eq. (1), is null-geodesically 
incomplete. Our proof also includes standard McVittie 
spacetimes as a particular case with constant rn{t), since 
the asymptotic behavior of geodesics is governed by the 
asymptotic properties of m{t) and H{t) through the func¬ 
tion ^(t), defined in Eq. (19), and its time derivatives. 

Therefore, in generalized McVittie spacetimes, the 
limit surface can be either a black hole event horizon, 
in part a black hole horizon and a white hole horizon, 
or even entirely a white hole horizon, depending again 
on the asymptotic properties of ^(f). We obtained this 
result by a similar method to the one used in [15], where 
the standard McVittie spacetime, with stricter assump¬ 
tions than those used in this paper, was proved to allow 
two of the three possible causal structures found here. 


since the white-hole-only case is not possible in standard 
McVittie spacetimes that obey our assumptions. Follow¬ 
ing the observation that the causal structures are distin¬ 
guished by the way ingoing null curves approach their 
limit at r = Too — either all from above r_(t), all from 
below r_ (t) or some from above and some from below — 
our proof relies on the fact that we only need to study 
the crossing of the geodesics with respect to a constant 
r = Too surface, instead of the time-dependent r_ (t) hori¬ 
zon. We then built two families of approximations for the 
ingoing null curves: one family that approximates from 
above and the other one from below such that we could 
write formally the solution to those approximations in 
order to determine whether there exists one curve that 
crosses the r = surface, or, alternatively, whether all 
of them fail to cross it. This led us to the condition on 
^(t) that determines the causal structure of almost all 
gMcVittie solutions (up to a set of zero measure) which 
satisfy the initial assumptions. 

We are provided with a large family of analytical solu¬ 
tions that are dynamical and can present evolving black 
holes and white holes isolated or in pairs. Those can 
be useful for many applications, as a laboratory to the 
study of the physics of dynamical black holes, toy mod¬ 
els for accreting astrophysical black holes, the study of 
bounded systems in cosmological backgrounds and, due 
to the richness of types of dynamical and causal horizons 
it presents, to the study of physics near different types 
of horizons, which are the scenario of many recent de¬ 
velopments such as the fluid-gravity correspondence in 
dynamical backgrounds [25], properties of vacuum solu¬ 
tions of modified gravity [12, 26], the search for the ther¬ 
modynamical laws of gravitational systems [27, 28] and 
even the debate on the proper definition of a black hole 
[29-31]. 
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Appendix A: On the approximation used for 

c(t) > C)(e-“"“‘) 

In section IIB, we claimed that the form of z{t) when 
^(t) > was given by (28), which we used to 

prove geodesic incompleteness in case 2. Here we show 
the deduction of this result. 

Integrating the first term of Eq. (22) by parts, we ob¬ 
tain 
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Jto Hoa L 


Hoa I 

c — aHot rt 


aH, 


0 Jto 


(Al) 


which we may write iteratively as 




n+1 


n=0 

-Coe 


Hoa J dt 


d” 


— aHni 


(A2) 


j 0 j 

where = / and 


Co = 

n—0 


( 1 


n+1 


\Hoa J dt 


3A7^(^o) ■ (A3) 


To prove this, we define y = and g{y) = Then, 

we have ^{t) = —y^g'{y)- The limit (B2) is written as 


lim 

y-fO 


y^g'jy) 

giy) 


(B3) 


We can express the denominator as 


giy) = y ^^y ^_Q^ = y ig'iy) + ^(y )). (B4) 


where limy_>.o ^(j/) = 0. Inserting Eq. (B4) into (B3), we 
obtain 


lim 

y-s-o 


yg'iy) 


y 

g'iy) + Uy) 

y-40 

1 + AM. 

^ + g'(y) 


(B5) 


which vanishes as j/ —)■ 0 unless limj,_j.o = — 1- This 
is the first result. 

In the special case 


Thus, we can write the flow of Eq. (21) as 
z(t)=(Zo-i?ooCo)e-“^“* 

OO / 1 \ "^ + 1 

+ i?ooE(-l)”^ 

n—0 

+ o((5). 


\Hoa 


d” 


(A4) 


Since Eq. (A2) is meaninful only when the series con¬ 
verges, we use this result in cases where the ^ domi¬ 
nates ^ and, consequently, ^ dominates higher deriva¬ 
tives of Thus, the leading term at late times is always 
the one with n = 0, and the first subleading term is 
n = 1. These assumptions corresponds to cases in which 
^(t) = o(e““^“0 (see Appendix B). The exceptions to 
this rule decrease either exponentially or faster and were 
treated separately in cases 3 and 1, respectively. 


^iy) 

lim 

g iy) 


= -1 


(B6) 


we cannot tell the value of the limit (B2) by the method 
above. We show here that one of three cases may happen: 


i. limt_j.oo 

ii. limt_j.oo 

iii. limt_j.oo 


id) 

id) 

id) 

id) 

id) 

id) 


0 ; 

L^O; 

OO. 


We have written = g'{y) + ^(y). Thus, ^(y) = 
^ - y'(y) and we have - 1. Therefore, 

Eq. (B6) is equivalent to 


Appendix B: When !^(t)| > as t —>■ oo 


giy) 


lim 

yg'iy) 


= 0 . 


(B7) 


In this Appendix we show why all the cases where ap¬ 
proximation (28) is valid are of slow decreasing, that is, 
correspond to case 2 in the analysis of section IIB. We 
assume ^ continuously differentiable Vf > 0, there exists 
T such that ^ 0, Wt > T and 


We can build functions with this asymptotic property by 
making 


g'iy) 


gjy) 
yhiy) ’ 


(B8) 


lime(f)=0. (Bl) 


Thus, limt_>.oo C(i) = 0- We will show that unless some 
specific conditions are met, we have 


lim 



(B2) 


where h{y) is any differentiable function such that 

lim h{y) = 0 . (B9) 

y-J-O 

This construction gives us y(y) of the form 

<,fa) = Cexp(|A). (BIO) 
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which implies 


^{t) = Cexp 



Y dt' 

i J 

t'h{^/t') 


iit) 

m 

m 




th{'^/t) 


dt' ' 


(Bll) 

(B12) 

(B13) 


Therefore, limt_>oo 


id) 

m 


is finite, unless asymptotically 


/i(Yt) < 0 ( 1 / 4 ). The first derivative ^ will only dominate 
^ in the case h{^/t) < 0(i/t), which implies, by (Bll), 
that ^(t) < 0(e“*)> which shows that those cases always 
fall in case 1. 


Appendix C: How determines the causal 
structure of gMcVittie spacetimes 


horizon we have d{t) < 0 and d{t) > 0, which means that 
it is possible for ingoing null geodesics coming from be¬ 
low to cross the r_ apparent horizon. If those curves do 
cross the horizon at a finite t, then they will reach the 
limit surface Too from below t-, since they cannot cross 
r_ back from above. Otherwise, if they do not cross r_ 
at finite t, they obviously reach the limit surface from 
above. From this reasoning, we can predict that, when 
f_ > 0 for large t, only the cases depicted in Figs. 4 and 
3a are possible. In other words, we need only focus our 
attention on geodesics approaching r_ from above. The 
distinction between these two remaining scenarios may 
be stated as the following: if there exists T such that, 
for any t > T, the ingoing geodesics just above r_ fail to 
cross the horizon, then our causal structure is like that 
in Fig. 3a. If there is no such T, in which case all ingoing 
geodesics coming from above the horizon do cross it in 
finite time, then our causal structure is that of Fig. 2. 

The next proposition is very useful to simplify the anal¬ 
ysis necessary to distinguish the remaining possible cases. 


In this Appendix, we show the steps that prove The¬ 
orem III.I. We have to consider separately the case 
f_ —>■ O’*' and f_{t) —>■ 0“. The proof is analogous for 
both cases, the only difference being in the sign of some 
functions. We then specialize to the case r_(t) > 0 for t 
sufficiently large. According to (42), this corresponds to 
^(t) —)■ 0+ as t —>• 00 . 

It is important to note that we need an additional as¬ 
sumption over ^(t) here, compared to Theorem II. 1. We 
assume that ^ (<) has a finite number of roots, such that 
it reaches zero either from above or from below and does 
not oscillate indefinitely between positive and negative 
values. More precisely: 

3T>0 such that ^(t) ^ 0, Wt > T, (Cl) 

which implies that there exists T such that either ^{t) >0 
and ^{t) < 0 for alH > T or ^(t) < 0 and ^{t) >0 for all 
t > T. 

Proposition C.l. If f-(t) > 0 for all t > T, and if 
r{t) is an ingoing null geodesics, governed by Eq. (10), 
that satisfy r(to) > r_(to) for some to > T, then we have 
r{t) > r-lf) for all t > to- 

Proof. We define d{t) = r^ft) — r{t) between the hori¬ 
zon and the geodesic r(t). Thus, for a given curve r{t), 
d(to) > 0. From the definition of r_, f(t) —0 as r —>■ r_, 
thus, for all t > to there exist e > 0 such that d{t) > 0 at 
d{t) = e > 0. Therefore, d{t) > 0 for all t > to- □ 

Proposition C.l implies that, if f_(t) > 0 for large 
times, ingoing null curves that are near r_ but below it 
remain so for indefinitely large times and reach r^o from 
below r_. This automatically excludes the causal struc¬ 
ture depicted in Fig. 2, as in that case there is no ingoing 
null curve that reaches the limit surface from below r_. 
However, for the curves that come from above, near the 


Proposition C.2. Let r(t) be an ingoing geodesic and 
r-{t) >0. If r satisfies r(t) > r-(t) for all t > to, then 

r(t) > Too, yt>to, (C2) 

and 

lim r{t) = roo • (C3) 

t—>oo 

Proof. We refer the reader to the proof of Proposition 
III-l of Ref. [15], but exchanging incresing sequences by 
decreasing ones and vice-versa. □ 

Thanks to Proposition C.2, one only has to analyze 

geodesics crossing roo, a fixed surface, eliminating the 
complication of having to consider the time-varying inner 
horizon r_(t). Moreover, it implies that the r^o surface is 
an accumulating point for geodesics coming from above, 
as we state in the following Corollary. 

Corollary C.3. If r(t) is an ingoing geodesic and 
r{to) > r_(to), then for all e > 0 such that the events 
[t-ifoo + e) liss in the regular region Vt > to, there exists 
t > to such that r(t) = r^o + £• 

This means that, given enough time, all geodesics ei¬ 
ther cross Too or reach values arbitrarily close to it. 
By Proposition C.2, every ingoing geodesic that never 
crosses r_ (t) never reaches r^o - Conversely, if an ingoing 
geodesic does traverse roo in a finite time interval, then it 
eventually crosses r_. Then, by studying only the neigh¬ 
borhood of Too, we may tell if geodesics do or do not cross 
the apparent horizon. With all this in mind, we can use 
the approximations made in Sec. IIB, which are valid for 
large t and r near roo, in order to use Eq. (18) and again 
defining z = r — roo, we obtain Eq. (21). Here we are 
interested in the error term of the approximation. We 
will use it to build approximations for the ingoing null 
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curves we are studying. We call them z±(t), which are 
curves that approximate z(t) from above and below re¬ 
spectively, valid when close to Too, that is, z < S. For 
an initial condition z(to) = zq > 0 and 0 < Xq < that 
is, above and close to the horizon, we are looking for two 
behaviors of the approximations z^(t) for large t: 

1. z'L{t) > 0 for all t > to. This guarantees that the 
ingoing geodesic does not cross the horizon in finite 
time. 

2. z'^{t) < 0 for some t > to. Then, its guaranteed 
that the ingoing geodesic does cross the horizon at 
a finite time t < t. 

Now, using the fact that there exists tr > 0 such that 
az > o(6) > —az , for zq < z < 0 , (C4) 

we build the approximations z'^{t) as solutions of 

= -aHoz%. + Roo^{t) ± az± , 

z‘^{to) = zo, (C5) 

such that we obtain, by Gronwall’s Lemma [32], 

z+{t) > zit) > z1{t), t>to, (C6) 

meaning that z'^{t) approximate z{t) from above and 
zZit) approximates z{t) from below as we intended. 

The solutions for z'^{t) are formally given by 

(C7) 

We remark that Zq is positive and the integral term 
is negative, since ^{t) < 0 for large t. Then we have to 
study the convergence of the functions F^{to,t) defined 


as 

Fl{to,t)= f < 0, (C8) 

Jto 

1. There exist L+{to) > 0 and ct > 0 such that 
limt^oo ^"+(^0,*) = L+{to) > 0; 

2. There exist cr > 0, such that limt_>oo FL(to, t) = 
-l-oo; 

which correspond to the following scenarios: 

Then, we conclude that 

1. If there exist tr > 0 such that F^{toR) —oo as 

t —>■ oo, then there exist an approximation from 
above such that the negative term of z'^{t) 

always can surpass the positive term and then z{t) 
changes sign. This implies that all curves z{t) 
changes sign and approach the limit surface from 
below and we obtain the structure depicted in fig. 4. 

2. If there exist cr > 0 such that FZ{to,t) converges 
to a value L_(to) < 0 as t —)■ oo, then, by a sim¬ 
ilar reasoning, there are approximations from be¬ 
low, (<) such that the positive term Zq is always 
larger in magnitude than the negative term, since 
we can make the integral term as small as we want 
by taking a larger to. Thus, in this case, we can 
say that there exists T > 0, such that for to > T, 
z{t) > 0 for all t > to. This correspond to the 
case where some curves approach the limit surface 
from below and some (the later ones) approach the 
limit surface from above. This correspond to the 
structure depicted in fig. 3a. 

This proves half of the Theorem III.l. The proof of 
the case f_ —)■ 0“ is identical but for the orientation of 
the z-axis. More details in this kind of argument can be 
found in Ref. [15]. 
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